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Abstract 

Ambiguities regarding the physics and the existence of the critical point 
(CP) on the QCD phase boundary still exist and the mist regarding the 
conjectured QCD phase boundary has not yet cleared. In this paper we 
extend our earlier study where we constructed a deconflning phase boundary 
using Gibbs' equilibrium conditions after using a quasiparticle equation of 
state (EOS) for quark gluon plasma (QGP) and an excluded volume EOS for 
the hadron gas (HG) and find the presence of a critical point on this phase 
boundary where the first order phase transition terminates. In this paper, 
we plot the difference in the normalized entropy density (s/T 3 ) between HG 
and QGP phases along the deconfining phase boundary and find that it 
vanishes at CP. Further we have shown the variation of the square of speed 
of sound (c^) for the HG and QGP separately and find that the difference 
(Ac^) between them along the deconfining phase boundary again vanishes 
at the CP of the boundary. We also plot the variation of shear viscosity 
to entropy density ratio {rj/s) in HG as well as in QGP phases separately 
with respect to temperature at different baryonic chemical potential (/is). 
The presence of a cusp like structure in r]/s again confirms the existence of 
CP on the deconfining phase boundary as predicted by some authors. These 
studies thus firmly indicate the precise location of CP as a second order phase 
transition point. 
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Precise mapping of the QCD phase boundary existing between hot, dense 
hadron gas (HG) and weakly interacting quark gluon plasma (QGP) and the 
location of a hypothesized QCD critical point (CP) have emerged as the 
most interesting and challenging problems before the experimental and the- 
oretical heavy- ion physicists today [1-3]. The possible existence of CP in 



Preprint submitted to Physics Letters B 



January 4, 2012 



the temperature (T) and baryon chemical potential plane of the QCD 
phase boundary was proposed a decade ago and it represents a second-order 
transition point where the first-order phase boundary terminates as T fur- 
ther increases and \ib decreases [3]. Its separation from the temperature 
axis (hb = 0) spans the region of a cross-over transition where mesons in 
this densely populated region gradually break-up due to the thermal fluctu- 
ations and are finally converted into a gas of quarks, anti-quarks and gluons 
without any kind of phase transition. Lattice QCD studies employing Monte 
Carlo simulation, fail at hb ^ because the absence of a probability measure 
precludes direct computations and hence we often use mathematical approx- 
imations in order to determine QCD phase diagram for non- vanishing values 
of hb- Some of these calculations have confirmed the existence of a cross-over 
chiral transition at fis = and T = 170 — 200 MeV and it ends at a CP with 
coordinates /x c /3T c < 1.0 [4-5]. However, other lattice calculations have cast 
a shadow of doubt on the existence of CP in the chiral phase transition [6]. 
Conclusively we still do not precisely know whether the conjectured phase 
boundary is an outcome of deconfinement and/or chiral symmetry restora- 
tion and whether the CP is related to the chiral phase transition alone [7]. 
Therefore, it seems interesting to study the nature of QCD phase transition 
and the existence of CP with the help of other phenomenological models. Re- 
cently we have constructed a first order deconfining phase boundary between 
HG and QGP phases employing Gibbs' criteria of thermodynamic equilib- 
rium in a hybrid model approach [8-9]. These investigations further revealed 
an interesting finding that the deconfining phase boundary between HG and 
QGP terminates at a CP beyond which solutions of Gibbs' conditions are 
not possible indicating the presence of a cross-over transition. The precise 
location of CP in (T, /ig) plane was obtained as (T c = 166, [i c = 155) when 
we used quasiparticle description for QGP [9] and this result is not much 
different from the location of CP obtained in the bag model EOS for QGP 
[8]. Our results get significant support from some earlier investigations where 
such hybrid model approaches [10-11] have been widely used. However, these 
studies do not reveal the precise nature e.g., the order of phase transition of 
CP. In this paper, we further extend our studies regarding the deconfining 
phase transition and the possible existence of CP by supportive evidences 
regarding the behaviour of the quantities such as normalized entropy den- 
sity s/T 3 , velocity of sound c s , and the ratio of the viscosity to the entropy 
density 77/s near the CP in our hybrid model. In particular we calculate the 
difference in entropy density As/T 3 as well as difference in sound velocity 
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Acj, between HG and QGP phases and we notice that these quantities vanish 
at the critical point and thus indicating a clear change in the order of the 
phase transition. By comparing the hadron-QGP transition with helium, ni- 
trogen and water at and near their phase transition points in the r]/s ratio, 
Csernai, Kapusta and McLerran [12] have recently shown that the variation 
of rj/s with temperature in both the phases can pinpoint the location of the 
critical point. The studies of Sasaki and Redlich [13] and Lacey et. al. [14] 
lend further support to this finding that rj/s ratio reveals a cusp like feature 
when we study its variation with T near the critical transition point between 
HG and QGP phases. Our aim in this paper is therefore, to explore further 
the precise location and the nature of CP on the deconfining phase boundary 
and to determine whether CP indeed is a point of second-order phase tran- 
sition. 

Recently we proposed a new thermodynamically consistent, excluded- 
volume model for the hot and dense HG [8-9, 15]. Our approach incorpo- 
rates the following new features. Besides thermodynamical consistency, our 
model uses full quantum statistics so that the phase boundary in the entire 
(T, hb) plane can be investigated without using any additional approxima- 
tion. Moreover, we incorporated excluded-volume correction arising due to 
a hard-core baryonic size only. We further assumed that the mesons can 
overlap and fuse into one another and hence cannot generate any hard-core 
repulsion. This is one big difference between our model and other models. 
We have demonstrated [15] that our calculations give a very good fit to the 
experimental ratios of multiplicities. Our model gives a proper description 
of all thermodynamical quantities including transport coefficients. The total 
pressure of the HG after incorporation of the excluded-volume effect can be 
written as [15-17]: 

phg =t(i-r)J2 hK + £ vT son - (!) 

i 3 

In Eq. (1), the second term on the right hand side gives the total pressure 
from all mesons in the HG having clearly a pointlike size. R = Yli nf x V® 
gives the fractional occupied volume due to all types of baryons and A« = 
exp(^r) is the fugacity of the particle. Jj represents the integral involving the 
distribution function : 

_ g_j r°° k A dk 1 
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where E { is the energy [Ei = ^Jk 2 + mf ) with /c as the momentum and mi 
as the mass of the ith baryon. 

Our calculation for the shear viscosity is completely based on the method 
of Gorenstein et al. [18]. According to molecular kinetic theory, we can write 
the dependence of the shear viscosity as follows [19]: 



n oc n I (|k|), 



(3) 



where n is the particle density, / is the mean free path, and hence the average 
thermal momentum of the baryons or antibaryons is: 



_ Jo k 2 dk k A 
j™k 2 dkA' 



(4) 



and A is the Fermi-Dirac distribution function for baryons (anti-baryons). 
For the mixture of particle species with different masses and with the same 
hard-core radius r, the shear viscosity can be calculated by using equation 
[18]: 

where n« is the number density of the ith species of baryons (or anti-baryons) 
and n is the total baryon density. 

To calculate the speed of sound at constant s/n, we have used the re- 
cent formulation of Cleymans and Worku [20]. The speed of sound at zero 
chemical potential is easy to calculate where it is sufficient to keep the tem- 
perature constant [21, 22]. However, the speed of sound (c s ) at finite chemical 
potential can be obtained by using the following extended expression [20]: 
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where the derivative dfj,B/dT and d\i s jdT can be evaluated by using two con- 
ditions, firstly by keeping s/n constant, and then imposing overall strangeness 
neutrality. Thus one gets [20]: 
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where L = nf + , is the sum of the strangeness density of baryons and 
mesons. Similarly R = nf + n™, stands for the sum of anti-strangeness den- 
sity of baryons and mesons. In all the above calculations, we have taken an 
equal volume V° = -^p- for each baryon with a hard-core radius r = 0.8 fm. 
We have taken all baryons and mesons and their resonances having masses 
upto 2GeV/c 2 in our calculation for the HG pressure. We have also used 
the condition of strangeness neutrality by putting ^ Si(nf — nf) = 0, where 
Si is the strangeness quantum number of the ith hadron, and nf(nf) is the 
strange (anti-strange) hadron density, respectively [23]. 

The EOS for QGP as used in this paper has been described in detail in 
Ref. [24, 25]. The detailed calcultaion regarding thermodynamical quanti- 
ties like pressure, energy density, particle density etc. can be found in our 
earlier work [9] . We have demonstrated that our model agrees well with the 
detailed features of the curves obtained in the lattice QCD calculations. We 
are thus confident in using our QCD quasi-particle even at finite [i B where 
lattice QCD fails miserably. We get a deconfining phase boundary after us- 
ing Maxwell's construction and we find that the phase boundary as well as 
CP obtained by this method matches well with those obtained in Bag model 
calculation [8]. Here we give the prescription used by us for calculating the 
transport properties e.g., shear viscosity, speed of sound etc. Our calculation 
for shear viscosity is based on the prescription used by Sasaki and Redlich 
[26] where they calculate the shear as well as bulk viscosity for QGP in the 
quasiparticle model. The shear viscosity in a medium composed of one type 
of particle/antiparticle can be obtained from the following expression [26] : 

1 f d 3 k k 4 

^ = 7^/(^)3^ k/o(l ± fo) + fff/o(l ± /o)] , (9) 

where k is the momentum and E = Vk 2 + M 2 is the energy with M being 
the thermal mass for the quark and gluon. ± is for fermion and boson (i.e., 
quark and gluon), respectively and fo (fo) is the equilibrium distribution 
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function for them, respectively and is given as : 

/o(/o) = (e (E ^ )/T ±l)-\ (10) 

where we use =F/x for f and fo, and ±1 for quark and gluon, respectively. In 
Eq. (9), r is the collision time which is determined by the thermal-averaged 
total scattering cross section (va) , with the relative velocity of two colliding 
particles v and the particle density n in equilibrium. Further g and g are 
the degeneracy factors for particle and antiparticle, respectively. We can get 
the shear viscosity for QGP simply by summing the Eq. (9) for all types of 
particles. The speed of sound for QGP is finally calculated by the following 
relation : 

where p is the pressure and e is the energy density of the QGP in our QPM. 

In Fig. 1, we demonstrate the results obtained for the trace anomaly 
factor (e — 3p) /T 4 in our hybrid model calculations using HG and QGP 
equation of state separately at hb = 0. We further compare our results with 
the results obtained in a recent lattice calculation [27]. We notice that our 
results yield an excellent fit to the lattice data. The success of our hybrid 
model which involves a separate and distinct description for both the phases 
(i.e., low temperature HG and large temperature QGP), is indeed excellent 
in reproducing the features of the lattice curve. Although such a combined 
description can still be treated as a crude one, but since lattice descriptions 
are still found quite unsatisfactory at 7^ 0, we emphasize that our hybrid 
model can provide an excellent phenomenological substitute for the formal 
theory. 

In Fig. 2, we have attempted to show what happens to the change in 
the entropy density at the CP of the phase diagram. We have calculated 
the difference ^§ = (s/T 3 )qqp — (s/T 3 )hg an d demonstrated its variation 
with respect to the coordinates of the phase transition points lying at the 
deconfining phase boundary. We find that ^§ ^ 0.0 and positive along 
the deconfining phase boundary in the case of first order phase transition 
which further indicates the presence of nonvanishing latent heat involved in 
the phase transition from HG to QGP. However, we surprisingly notice that 
^§ at the CP and thus the CP can be taken as a point where the first 
order phase boundary terminates and phase transition changes its order. 
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Figure 1: Variation of trace anomaly (e — 3p) /T 4 with respect to temperature. Solid line 
represents the low temperature behaviour of trace anomaly calculated in our excluded 
volume model for HG. Dash-dotted curve is the result obtained from quasiparticle model 
of QGP. Lattice data points are taken from Ref. [27]. 
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Figure 2: Variation of (As/T 3 ) = (s/T 3 )qgp — (s/T 3 )hg with respect to coordinates of 
various phase transition points on the (T,/jb) phase boundary. We have used transition 
points from Ref. [9]. 
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Figure 3: Variation of c 2 s )qgp (dash-dotted line ), (c 2 s )hg (dashed line), and (Ac 2 ) = 
(cDqgp ~ [c 2 s )hg (short-dashed line) with respect to coordinates of various phase tran- 
sition points on the (T,hb) phase boundary. We have used transition points from Ref. 
[9]. 
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Figure 4: (A) Variations of shear viscosity to entropy density ratio (r//s) for gluon plasma 
with respect to T/T c . Solid line is the result from our calculation. Long-dashed line is 
the result obtained in Ref. [35] and short-dashed line is taken from Ref. [36]. The lattice 
data with 16 3 x 8 and 24 3 x 8 lattice are from Rcfs. [37] (triangles and squares) and [38] 
(filled circles). (B) Variations of i]/s obtained in our calculations for HG and QGP with 
respect to temperature at various \i B . 
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This is certainly an interesting finding. Although we have used a hybrid 
model description for QGP and HG phases, the vanishing of net entropy 
density at the CP indicates that there exists a continuity in these two model 
descriptions. This means that the parameters in both descriptions had been 
suitably adjusted. 

Sound velocity is an important transport property of the matter created 
in nucleus-nucleus collision experiemnts because the hydrodynamic evolution 
of this matter strongly dependes on it. Speed of sound is related to the speed 
of small perturbations produced in its local rest frame. A minimum in the 
c s has also been interpreted in terms of a phase transition [20-21, 28-33] 
a large number of degrees of freedom is present which causes difficulty in 
the propogation of the sound wave in the medium. Further, Chojnacki and 
Florkowski [34] proposed that a shallow minimum in the speed of sound near 
the smooth joint of HG and QGP phases in a hybrid model corresponds to 
the presence of a cross-over transition. In Fig. 3, we have separately shown 
the variation of square of speed of sound i.e., c 2 for HG and also for QGP. We 
have also shown the difference (Ac;?) = (c 2 )qgp — (c 2 )hg an d demonstrated 
its variations with respect to the coordinates of the phase transition points 
lying at the deconfining phase boundary. We again find that Ac 2 . ~ 0.0 at the 
CP. Thus this results lends further support to our previous finding regarding 
the change of the order of the transition at the critical point. 

In Fig. 4 (A), we plot the variation of rj/s of gluon plasma with respect 
to T/T c , where T c is the critical temperature. We compare our model result 
with the results obtained by Khvorostukin et. al. [35] and Bluhm et. al. 
[36]. We have also shown a comparison with the lattice gauge calculations 
[37, 38]. We observe that the results obtained in our calculation agrees well 
with the lattice data even at large temperatures. However, the error in lattice 
simulation result as given in Fig. 4 (A) is quiet large. Our result lies in close 
aggreement with the result obtained by Khvorostukin et. al [35]. This result 
naturally gives us additional confidence to calculate rj/s of QGP by using 
QPM. In Fig. 4 (B), we plot the variation of shear viscosity to entropy 
density ratio for HG and for QGP, respectively with respect to temperature 
at fiB — 300 MeV (dashed-line) and at critical potential \x c = 155 MeV (solid- 
line). At /ib — 300 MeV, we observe a discontinuity in rj/s at the joint of both 
the curves for phases i.e., HG and QGP, at first order phase transition point. 
Further we observe an upward jump in rj/s going from low temperature HG 
phase to high temperature QGP phase which is in accordance with the result 
obtained by Sasaki and Redlich [13]. At critical chemical potential /i c = 155 
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MeV and temperature T c = 166 MeV [9] we get a cusp like behaviour in rj/s 
while going from HG to QGP phase as shown in Ref. [13, 39]. 

Thus above results give a firm indication that the order of phase tran- 
sition changes at CP for the deconfinement phase transition. It should be 
added here that many authors in the past have used two separate equations 
of state for QGP and HG and obtained a tentative explanation for an an- 
alytic and smooth cross-over and CP in their models [10-11]. Our model 
presents a similar picture. The physical mechanism involved in this calcu- 
lation is intuitively analogous to the percolation model where a first order 
phase transition results due to 'jamming' of baryons which thus restricts the 
mobility of baryons [40-41]. However, in the percolation model we do not 
have any comparison to what one gets in the QGP picture. Here we ex- 
plicitly and separately consider both the phases, i.e., HG as well as QGP 
and hence it gives a clear understanding how a first-order deconfining phase 
transition can be constructed in nature. At low baryon density, overlapping 
mesons fuse into each other and form a large bag or cluster, whereas at high 
baryon density, hard-core repulsion among baryons, restricts the mobility of 
baryons. Consequently we consider two distinct limiting regimes of HG, one 
is a meson-dominant regime and the other is a baryon dominant region and 
the point joining the two is CP in our models. 

In conclusion, searching for the precise location and the nature of the crit- 
ical point (CP) on the QCD phase diagram still poses a challenging problem. 
Although various calculations have predicted its existence but the quanti- 
tative predictions regarding its location wildly differ. Experiments face an 
uphill task in probing the CP in QCD phase diagram because a clarity in 
theoretical prediction is missing. Moreover, many unstudied problems such 
as short lifetime and the reduced volume of the QGP formed at colliders also 
affect the location of CP and its verification [42]. In these circumstances, 
we believe that our results will throw light on the nature of CP existing in 
deconfining phase transition. 
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